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We present percolation strategies based on multipartite measurements to propagate entanglement 
in quantum networks. We consider networks spanned on regular lattices whose bonds correspond to 
pure but non-maximally entangled pairs of qubits, with any quantum operation allowed at the nodes. 
Despite significant effort in the past, improvements over naive (classical) percolation strategies have 
been found for only few lattices, often with restrictions on the initial amount of entanglement in 
the bonds. In contrast, multipartite entanglement percolation outperform the classical percolation 
protocols, as well as all previously known quantum ones, over the entire range of initial entanglement 
and for every lattice that we considered. 

PACS numbers: 03.67.Bg, 64.60.ah 



Quantum networks, which consist of neighboring nodes 
(or stations) sharing entangled pairs of particles, are the 
future of quantum communication [l[. Establishment of 
entanglement between two nodes of the network that are 
separated by an arbitrarily large distance represents the 
most important task and challenge in the study of quan- 
tum networks. The first solution to this problem was by 
means of quantum repeaters [2] , a (one-dimensional) line 
of intermediate stations that share multiple noisy entan- 
gled pairs. Long-distance entanglement is then achieved 
via entanglement distillation at the repeater stations and 
the use of quantum memories. 

A new approach for long-distance entanglement distri- 
bution was introduced in [3J, exploiting the connectivity 
of multi-dimensional networks. In this setting networks 
are spanned on regular lattices, where bonds are par- 
tially entangled pure states of two qubits and with arbi- 
trary quantum operations allowed at the nodes. Using lo- 
cal operations and classical communication (LOCC) one 
can probabilistically convert such states into perfect sin- 
glets Therefore, the problem is solved by classical 
bond percolation: entanglement propagation is possible 
when the conversion probability exceeds the percolation 
threshold of the lattice 0, 0|. However, it has been 
shown in @, [1| that this classical entanglement per- 
colation (CEP) is not optimal. In fact it is sometimes 
advantageous to consider quantum entanglement perco- 
lation (QEP): one precedes the singlet conversion by a 
suitably designed set of local quantum operations, and 
applies CEP afterwards. Despite several studies devoted 
to these strategies, little is known about entanglement 
propagation protocols in general lattices. 

In this Letter we introduce a powerful class of QEP 
protocols that exploit multipartite entanglement. In con- 
trast to the strategies studied so far, which solely employ 
entanglement swapping and conversion into singlets, we 
make full use of both the classical and quantum aspects 
of quantum networks: connectivity and multipartite en- 



tanglement, respectively. The interplay of geometrical 
lattice transformations and entanglement manipulations 
is in fact a key ingredient to surpass CEP. Before ex- 
plaining our protocols in detail, let us first state the main 
results of our work: i) multipartite strategies relate for 
the first time entanglement propagation to classical site 
percolation in a natural way [9|; ii) they systematically 
outperform all previously known classical or quantum 
strategies, regardless of the initial entanglement of the 
bonds and for every lattice that we considered; iii) they 
are obviously applicable to any lattice relevant for phys- 
ical problems. 

Let us now be more specific and recall the main idea 
of CEP. Each bond of the lattice is given by the state 
\p) = VPo |00) + V /^T |H), with (p > ipx and <p +pi = 1, 
and can be optimally converted into a singlet (or equiva- 
lently into the Bell pair ]$+) cx |00) + |11», with prob- 
ability p(ip) = 2tpi. If p(ip) is larger than the thresh- 
old for the corresponding bond percolation, then a giant 
cluster of maximally entangled pairs appears and long- 
distance entanglement between two nodes of this cluster 
is achieved by performing entanglement swappings. 

Concerning multipartite protocols, the main idea is 
to produce Greenberger-Horne-Zeilinger (GHZ) states, 
which are the generalization of the Bell pair |<E> + ) to n 
qubits: |GHZ„) cx |00 . . . 0) + 1 11 ... 1) . For instance, 
let us show how to build a tripartite GHZ state from 
two neighboring bonds, that is, from two copies of \<p) 
sharing a common node. Guided by the optimality of 
Bell measurements in the ZZ basis [id ] for one-repeater 
entanglement swapping @, Sect. Ill A], we perform on 
the two qubits of the common node a measurement M 
consisting of the two operators |0)(01| + |1)(10| and 
|0)(00| + |1)(11|. The first outcome leads to a perfect 
GHZ state among the three nodes, but with probability 
Po + Pi t ne second outcome yields the (unnormalized) 
state <po |000) + ipi |1H)- As for a pair of qubits, the 
latter state is transformed into a GHZ state with prob- 
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FIG. 1: (Color online) From bond percolation to site percolation on a lattice with higher connectivity. First Step: Each 
measurement probabilistically transforms two links into a GHZ state on three nodes, which we depict by a triangle in £' and 
a site in C. Second step: two neighboring triangles can be merged into a bigger connected structure, which we represent by an 
edge in the transformed lattice, a) Measurement pattern to transform the (4,8 2 ) lattice into a lattice with a better threshold. 
The right panel shows the corresponding site percolation problem; b) The hexagonal lattice is transformed into a lattice that can 
be denoted by 1 /3(A-, 5 2 ) + 1 /3{3, 4, 5, 4) + 1 /3(3, 4, 3, 5 2 ); c) Measurement patterns for the dice and (3 2 , 4, 3, 4) lattices, resulting 
in non-planar graphs (in the latter construction some links are converted into Bell pairs as in CEP). Semi-transparent circles 
highlight the nodes with highest coordination number used for computing 9. 



ability 2(p±/{ip^ + <p\). Summing the two possibilities, 
we find that two copies of \<p) are converted into a GHZ 
state of three qubits with optimal probability 2tp±. The 
multipartite method derives its power from the fact that, 
while the value of this probability is the same as for the 
two-qubit swapping used in previous entanglement perco- 
lation studies, we now have three entangled qubits rather 
than two. 

Although in this work we use only three-qubit GHZ 
states as building blocks for multipartite entanglement 
percolation, we can generalize the above procedure to 
construct GHZ states of n + 1 qubits starting from n 
copies of \tp) sharing a common node. In this case 
we have 2™ _1 measurement operators E m of the form 
|0)(m| + |l)(m| , where fh is the complement of m writ- 
ten in base 2. One can check that an (n + l)-qubit GHZ 
state is then created with probability 



P(GHZ 
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where [x] denotes the integer part of x. This probability, 
which is only optimal for n < 2, can be understood as 
a sequence of n — 1 entanglement swappings in the ZZ 
basis on a chain of n links, as calculated in [7|, Eq. D3]. 

The multipartite strategy consists in creating a lattice 
C, where nodes represent the GHZ states created from 
a lattice CJ Two vertices in C are connected by a 



bond if the corresponding GHZ states share a common 
node in the original lattice, see Fig.QJi,. This defines a site 
percolation process with occupation probability p = 2<pi. 
Above the percolation threshold of the new lattice, en- 
tanglement is propagated over a large distance as follows. 
Consider the situation in which two GHZ states of size n 
and m sharing one node have been created. One builds a 
larger GHZ state on n+m— 1 particles with unit probabil- 
ity by applying the measurement M on the two qubits of 
the common node. This generalized entanglement swap- 
ping can be iterated, and one eventually gets a giant GHZ 



state spanning the network. Finally, note that, given a 
GHZ state of any size, a perfect Bell pair is created be- 
tween any two of its qubits by measuring all other qubits 
in the X basis. 

At this point let us define the quantities characterizing 
long-distance entanglement. Connectivity in percolation 
is characterized by the correlation length £(p) which ap- 
proaches as p — ► and as p — ► 1 , and diverges as p j p c 
and as p f p c . In CEP we are interested in the prob- 
ability P'(A «-> B) of creating a Bell pair between two 
nodes A and B separated by a distance L and chosen 
from the set of all nodes with the largest coordination 
number Z max - We shall assume that A and B satisfy 
L » £'(p)- For p < p' c , P'(A B) decays exponen- 
tially in L/£'(p), while for p > p' c the two nodes are con- 
nected only if they are both in C In this limit the events 
{A £ C'} and {B E C'} are independent, which, together 
with translational invariance, gives P'(A *-> B) = 0' 2 (p), 
where 0'(p) = P'{A € C). Thus the problem is reduced 
to studying 9'(p). For QEP, we further restrict the two 
distant nodes to be chosen only from those at which no 
measurement M is to be made. In fact only these nodes 
retain the coordination number Z max . In terms of the 
probability measure on the site percolation process on C 
we have 9{p) = P(Ui{Ai 6 C}), where the union is over 
all the Z max sites A4 possessing a qubit that is also in 
A. It follows that P(A <-> B) = 9 2 {p), and thus a direct 
comparison of 9'(p) and 9(p) tells us which from CEP or 
QEP is favorable. Note that this is the case because in 
all our examples we have p = p' = p = 2(p\. 

Now that the basic ingredients have been presented, we 
apply them to several natural examples of lattices. In the 
following paragraphs we prove for eight lattices that mul- 
tipartite entanglement percolation protocols beats CEP. 
First, we examine the thresholds and show that p c < p' c 
for each lattice. Then, we compute the expansions 9(p) 
and 9'(p) in the high-density (i.e. maximally- entangled) 
limit to prove that 9(p) > 9'(p) as p tends to unity. 
Finally, using Monte Carlo techniques, we show that 
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TABLE I: Improvement of entanglement percolation thresh- 
olds using a multipartite strategy, and relative gain A = 
1 — Pc/p'c- We performed Monte Carlo simulations to cal- 
culate pc for lattices 2-8; all other values can be found, with 
higher precision, in [1, 0, 0- 



P(A <->£?)> P'(A <-*• B) for all p > p c . 

We start by considering the Archimedean (4, 8 2 ) lat- 
tice, see for instance [3] for the notation. In Fig. 
we propose a measurement pattern indicating on which 
qubits the measurement M is applied: the new lattice 
is denoted by 2 / 3 {3 2 , 6f) + 73(3,6,3,6) and its critical 
point is p c « 0.6499 [12j. Since the original threshold 
is p' c w 0.6768 the proposed strategy yields an im- 
provement over CEP. It is interesting to note that the 
transformed lattice, arising from simple quantum oper- 
ations, is among the rather exotic examples considered 
in other studies. For example, this two-uniform lattice is 
121 ] , where a quantitative relation between 
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percolation thresholds and the Euler characteristic is put 
in evidence. 

We studied many other lattices and found that multi- 
partite entanglement percolation leads to better entan- 
glement thresholds than CEP in every case, which sug- 
gests that this may be a universal result. To support this 
fact, let us describe seven more well-known lattices and 
the transformations that improve their critical points. In 
Fig. [TJ3-C we show the measurement patterns for three of 
these lattices, namely the hexagonal, dice and (3 2 , 4, 3, 4) 
lattices. The four other constructions are similar to the 
latter one: first, one or two overlapping square lattices of 
double size are created from the original lattice. Then, 
the remaining links are paired in a regular way, as de- 
scribed in @, Fig. 14] and @, Fig. 7]. The key point of 
our constructions, with respect to previous protocols, is 
that we do not get disjoint lattices anymore, but rather 
connected ones since the middle qubits can be used to 
propagate entanglement through the network. We did 
not find published critical values for site percolation in 
the resulting lattices, mainly because they are non-planar 
or non-regular graphs. Therefore, we turned to Monte 
Carlo simulations to provide numerical results. The val- 
ues obtained are summarized in Tab. HI thresholds are 
not only better for all lattices, but the gain is often sig- 
nificant, especially for lattices of high connectivity. 



TABLE II: Series expansions of 0'(p) and 0(p), for p = 1 — 
e and < e -C 1. These formulae have been derived for 
a fraction / = d/d! of nodes, where d! (d) is the density 
of nodes of higher connectivity in the original (transformed) 
lattice. Note that d! = 1 for the Archimedean lattices, since 
by definition all their vertices are equivalent, while d! — 1/3 
for the dice and d! — 1/2 for the bowtie lattice. 



We provided numerical evidence that multipartite en- 
tanglement percolation yields better thresholds than 
CEP, and therefore that the connection probability be- 
tween two widely separated nodes is increased when the 
bond entanglement p(ip) is small and lies in the interval 
(PaPc)- We now consider the opposite regime in which 
connections are highly entangled, with p = 1 — e and 
< e <C 1 . In this situation an analytical study becomes 
possible, and we use the perimeter method [8j, App. C] 
to compute high-density series expansions of 9'{p) and 
9{p). The two lowest non-trivial orders of the expansions 
are easily calculated and are sufficient to show that our 
strategy leads to an improvement over CEP for all lat- 
tices, see Tab. [TTJ Interestingly, this improvement does 
not appear at the first non-trivial order since there is no 
way to increase the number of independent connections 
of a node. In fact generalized entanglement swappings 
remove such connections, so that only the nodes where 
no quantum operation is performed share the property 
6{p) > 0> (p). 

The previous analyses show that the proposed QEP 
strategy works well near the classical thresholds and near 
the ideal situation of perfect connections. We thus expect 
that our transformations give P(A <->£?) > P'(A B) 
for all values of bond entanglement. We performed Monte 
Carlo studies for the eight lattices studied in this work, 
using importance sampling via the perimeter method 
in the high-density region, attesting that 9{p) is indeed 
strictly larger than 9'(p) for p £ (p c ,l). This is shown 
for three lattices in Fig. [51 

It is now legitimate to wonder about the optimality of 
our multipartite entanglement schemes. It is not diffi- 
cult to see that they cannot be optimal for every value of 
the bond entanglement. In fact, let us consider a square 
lattice and the deterministic "centipede" strategy intro- 
duced in 

0, Sec. VI, A]. The idea is that part of the 
links are used to create a spiral centipede whose "legs" 
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FIG. 2: Numerical determination of 9 — 9' showing that mul- 
tipartite percolation yields better long-distance entanglement 
for any initial entangled state \cp), with p = 2<pi. We show 
here the results for the triangle, square and hexagonal lat- 
tices (from left to right). High-density expansions are shown 
as light grey lines. 



have length one. To that end, one performs two entan- 
glement swappings (measurements in the XZ basis) on 
three links of a square. The resulting state is combined 
with the fourth unused link to get a maximally entan- 
gled state [4J. Clearly, this method can only be applied 
to lattices whose nodes have four or more neighbors. For 
square lattices, it leads to / = 1/2 and 9{p) = 1 above 
the point p c « 0.684. Therefore, in the case of highly 
entangled bonds, this methods is better than the present 
protocol, which gives 9(p) = 1 only at the trivial point 
p — 1. However the centipede strategy fails whenever 
p lies between 0.5 and p c , and therefore does not even 
reach the classical threshold. 

In conclusion, we have proposed a new approach to en- 
tanglement percolation in quantum networks that makes 
significant improvements over CEP. It makes use of mul- 
tipartite entangled states to improve the probability of 
long-distance entanglement generation between nodes, 
regardless of the lattice topology or the degree of en- 
tanglement of the bonds. In fact, we have provided both 
analytical and numerical evidence for several well-studied 
lattices showing that our protocols yield better results 
(lower thresholds and higher probabilities) than all pre- 
viously known strategies. At the moment, we do not 
know if the thresholds obtained from our multipartite 
QEP arc optimal. This brings us to the important open 
question of whether there is a minimum amount of ini- 
tial pure-state entanglement necessary for long-distance 
entanglement. In one-dimensional chains we know that 
this is possible only if the connections are maximally 



entangled, but while any percolation protocol defines a 
sufficient condition for long-distance entanglement dis- 
tribution in higher-dimensional systems, a proof of the 
necessary amount of entanglement is lacking. 

Finally, the most important challenge is to generalize 
the results obtained so far to the realistic case of mixed 
states. Recently, some results have been presented in [lij ]. 
where bonds are mixed states of rank two, which may 
be a first step in this direction. However realistic noisy 
networks are best described by full-rank mixed states, 
and even if we do not touch this problem in this Letter, 
the generality of our results suggests that multipartite 
strategies may also be of great importance in that case. 
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